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THE FFT: 


Fundamentals and Concepts 


With the continuing trend toward more economical computer implementation of 
digital techniques, we are beginning to see new ways of looking at old things. The 
Fourier transform, a mathematical tool for expressing time domain functions in the 
frequency domain, has been available to us since the early 1800's; but the complex- 
ity of applying it to anything other than the simpler functions has, until recently, 
limited its usefulness. The fast Fourier transform, (FFT), a digital version of the 
Fourier transform, has eased transformation between the time and frequency domains 
and provided the vehicle for carrying Fourier techniques to Many new areas of 
application. 

Waveform sampling and digitizing, along with the speed and accuracy of computer 
techniques, have simplified the task of waveform processing in either the time or 
frequency domains. With the TEKTRONIX® Digital Processing Oscilloscope, fast 
Fourier transformation of a waveform is reduced to a matter of obtaining an oscil- 
loscope representation of the waveform and pushing a few buttons. The complexities 
of obtaining the transform are delegated to a minicomputer and some powerful soft- 
ware; however the responsibilities of understanding the process and interpreting the 
results still rest with us, the user. To this end, understanding the process and 
interpreting the results, let us direct our attention first to the time and fre- 
quency domains and then to the Fourier transform and some of the approximations 


necessary to its implementation through the fast Fourier transform. 





The Time And Frequency Domains 


To begin with, let us look at a waveform in the time domain and see how it 
relates to the frequency domain. Then we can look at it in the frequency domain 


and see how it relates back to the time domain. 
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Figure 1 


Figure 1 shows three different representations of the same waveform. In l-a, 
we see a periodic nonsinusoidal waveform in the time domain. In l-b, we see another 
time domain representation of the same waveform. This representation is not of the 
waveform itself, but of the amplitude-phase relationship of the waveform's frequency 
components. Summing the amplitudes of each frequency component in l-b would bring 
us back to the representation of l-a. The waveform's frequency domain representation 
is shown in l-c. 

Let us go back to l-a again. As was mentioned before, this is a periodic non- 
sinusoidal waveform in the time domain.. One interesting and useful property of 
periodic waveforms is that they are all composed of sine waves. That is, any peri- 
odic nonsinusoidal waveform may be generated by adding together various sine waves 
of appropriate frequency, amplitude, and phase. Figure l-b contains the sine waves 
appropriate to generating the waveform of l-a. Figure l-c, sometimes referred to as 
a line spectrum, contains nearly the same information as l-b. The presence of each 
frequency component is indicated by a vertical line. This vertical line's location 
on the abscissa indicates the frequency of the component, and its magnitude is 


indicated by the height of the line. The line indicates that energy exists only at 
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the frequency where the line is located. Phase information could be contained on 
another diagram; however, for the case at hand, each component has zero phase shift 
with respect to the others so a frequency domain diagram of phase would be rather 
unexciting. 

Now the question arises, "How do you get from a to b?" In other words, we 


have only the waveform shown in Figure l-a, how do we find its frequency components? 


The Fourier Transform 


The fundamental frequency, indicated as f, in Figure l-c and the lowest fre- 


quency sine wave in Figure l-b, can be ау very easily by measuring the period 
of the waveform in Figure l-a and taking its reciprocal. But, what about the rest 
of the frequency components? What is their frequency, their relative magnitude, and 
their phase? These same questions bothered a French mathematician, J. B. J. Fourier 
(1758-1830), and as a result we have today a transform bearing his name. 

The Fourier transform is simply a mathematical tool that allows us to take a 
function of time, f(t), and re-express it as a function of frequency. When we 
transform a continuous periodic function of time, such as the waveform shown in 
Figure l-a, we will obtain a series of weighted impulse functions in the frequency 


domain, such as shown in Figure l-c. The weight assigned to each impulse function 


will correspond to the amplitude of that frequency component. To obtain a basic 


understanding of what has just been said, let us consider the case shown in Figure 2. 
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Figure 2-a shows a periodic waveform, a sine wave, that is continuous. By 
continuous, we mean that the function exists at all points in time from minus infin- 
ity to plus infinity. If we Fourier transform the continuous sine wave of Figure 
2-a to the frequency domain, we will see the frequency domain representation of the 
sine wave as the weighted impulses in Figure 2-b. A weighted impulse is simply a 
"spike" of zero width that has had its amplitude adjusted, or weighted, to reflect 
some quantity. In the case of Figure 2-b, the impulses are depicted by lines whose 
amplitude reflects the peak amplitude of the sine wave that was Fourier transformed. 
The positioning of these weighted impulses on the frequency axis serves to indicate 
the frequency of the sine wave. In Figure 2-b, note that the single frequency of 
the sine wave has been indicated twice, once in the positive frequency domain and 
once in the negative frequency domain. This does not mean that the sine wave is 
composed of two frequencies. Negative frequencies are just a mathematical dupli- 
cation or mirror imaging of the positive frequency domain in the negative frequency 
domain. In some cases, the negative frequency domain is necessary to complete wave- 
form description; however, for the most part, we can obtain whatever information we 
need from the positive frequency domain. 

Except for serving the function of defining basic terms and conveying the 
concept of transforming between time and frequency, the case shown in Figure 2 is 
really trivial. After all, we can easily measure the peak amplitude of a sine 
wave and the reciprocal of the period will give us frequency. We do not have to 
worry about harmonics because a sine wave does not contain harmonics (harmonics are 
frequency components that are multiples of the fundamental). A not so trivial case, 


a square wave, is illustrated in Figure 3. 
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Figure 3-a shows a square wave that we will specify as being continuous for 
all time; that is, it is defined for all points in time and is not time limited. 
Because it is periodic its spectrum will consist of weighted impulses as shown in 
Figure 3-b. The spectrum of Figure 3-b is known as a "discrete-frequency spectrum" 
or "line spectrum" and will be the subject of later discussion. As was previously 


mentioned, negative frequency components are simply duplicates of the positive fre- 


quency components so we will confine our discussion to the positive frequency domain. . 


In Figure 3-b, the weighted impulse occurring just to the right of zero fre- 
quency, the one having maximum magnitude, is referred to as the fundamental. Its 
amplitude corresponds to the peak value of the square wave and its frequency cor- 
responds to the reciprocal of the square wave's period. The next weighted impulse 
in Figure 3-b, occurring just to the right of the fundamental, has an amplitude 
less than that of the fundamental and a frequency greater than that of the funda- 
mental. Іп fact, theory shows that for the square wave of Figure 3-a the amplitude 
of the second impulse will be one third of the fundamental amplitude and the fre- 
quency will be three times that of the fundamental. Continuing with theory, the 


‘amplitude of the next impulse will be one fifth that of the fundamental and its 


frequency will be five times the fundamental. The next impulse will be related to 


the fundamental in the same manner by one seventh amplitude and seven times the 


‘frequency and so on until we have an infinite number of frequency terms. 


All of the weighted impulses to the right of the fundamental occur at integer 
multiples of the fundamental frequency and are referred to as harmonics. The fact 
that the harmonics in Figure 3-b occur at odd integer multiples of the fundamental-- 
3, 5, 7, ...--gives rise to their being referred to as odd harmonics. Had the har- 
monics occurred at even integer multiples of the fundamental--2, 4, 6, ...--they 
would be referred to as even harmonics. 

So far, we have seen that the Fourier transform may be used to describe a 


periodic waveform in terms of its frequency components and their magnitude. With 


frequency and magnitude determined, the phase relationship of the frequency compo- 


nents is the only remaining information required to completely specify a waveform 
in the frequency domain. The Fourier transform does provide phase informàtion, and 
if we look at Figure 4, we can see the phase spectrum that would be obtained for 


the square wave of Figure 3-a. 
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Figure 4. The phase spectrum obtained by 
Fourier Transformation of the 
square wave in Figure 3-a. 


Figure 4 contains lines on the frequency axis corresponding to the spectral 
lines in Figure 3-a. The magnitude of each line in Figure 4 serves to indicate the 
phase of its associated frequency component. Inspection of the positive frequencies 
in Figure 4 shows that each harmonic of the square wave has been assigned a phase of 
7/2 (90 degrees). Since each positive component has the same phase angle, it follows 
that there is zero phase shift between the components. These observations, with a 
change of sign, are also applicable to the negative frequency spectrum. In comparing 
the positive and negative frequency phase spectra, we can see that there is T (180 
degrees) phase shift between the positive frequency components and the negative 
frequency components. With frequency domain magnitude and phase information, like 
that shown in Figures 3-b and 4 respectively but extending in frequency from nega- 
tive to positive infinity, the time domain square wave of Figure 3-a is completely 
specified in the frequency domain. | 

By now we should realize that we can use the Fourier transform to get from the 
time domain to the frequency domain. We can now ask, "How do we get back to where 
we started?" The answer is simple - use the inverse Fourier transform. 

If a waveform is completely specified in the frequency domain, magnitude and 
phase provided, application of the inverse Fourier transform will re-express the 
waveform in the time domain. In other words, inverse Fourier transformation of the 
information in Figure 4 and Figure 3-b will result in the waveform of Figure 3-a. 

Up until now, the Fourier transform has been discussed in terms of magnitude 
and phase. The point should be made here that the Fourier transform can also be 


expressed in a complex form. That is, rather than specifying magnitude and phase, 
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ме сап specify real and imaginary components in the frequency domain. This is not 
too difficult to understand when we look at the analogous case of an RC impedance 
where we can specify the impedance as a magnitude and an angle (26/9) ог аз а геа1 
and imaginary part (Rt5X,)- 

Let us continue our discussion of the Fourier transform by looking at a repet- 
itive square pulse and observing the effect of increasing its period without varying 


its duration. The repetitive square pulse for consideration is shown in Figure 5-a. 
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magnitude of the repetitive 
square pulse in a. 


Figure 5 


As a first step, compare Figure 5-a with Figure 3-a. In so doing, we can 
quickly see that, except for the zero reference, the repetitive square pulse of 
Figure 5-a is identical to the square wave in Figure 3-a. The only difference 
between the two is that the repetitive pulse has no negative going portion. 

Next, let us make the stipulation, as we have in past discussion of the Fourier 
transform, that the waveform in Figure 5-a be continuous as shown from time equal 
minus infinity to time equal plus infinity. With this stipulation fulfilled, the 
Fourier transform magnitude of the repetitive square pulse will appear as the line 
spectrum in Figure 5-b. 

Since we previously observed that the repetitive square pulse in Figure 5-a was 
quite similar to the square wave in Figure 3-a, we might expect to see the simi- 
larity carried over into the frequency domain magnitude. In other words the fre- 
quency content might reasonably be expected to be the same; however we should not 
be too optimistic about similarity in phase. If we compare Figures 5-b and 3-a, 


the Fourier transform magnitude for the repetitive square pulse and the square wave 
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respectively, we will see that both have essentially the same frequency components, 
although at different magnitudes. The major difference between the two occurs at 
zero frequency, dc, where we see that the repetitive square pulse makes a major 
contribution that does not exist for the square wave. 

We should not be too surprised at the dc component appearing in Figure 5-Ь. 
The waveform that was transformed consisted of pulsating dc, in fact, if we look 
closely we can see that it was at full value half of the time and at zero the rest 
of the time (50 per-cent duty cycle). By inspection, we can see that the average dc 
value of the repetitive square pulse will be equal to half of the pulse amplitude. 
If we look at the Fourier transform magnitude in Figure 5-b, we can see that the dc 
component in the frequency domain corresponds to the average dc in the time domain 
waveform. | 

Let us increase the period of the pulse train shown in Figure 5-a until we 
obtain the pulse train shown in Figure 6-a. The effect of this increased period on 


the frequency domain is shown in Figure 6-b. 
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a.) A repetitive square pulse. b.) The Fourier transform 
magnitude of the repetitive 
square pulse in a. 


Figure 6 


The pulse train or repetitive square pulse in Figure 6-a still has an average 
dc value, but since the period has become longer than the pulse duration we will see 
a lower average dc than was seen in Figure 5-a. The Fourier transform magnitude, 
shown in Figure 6-b, has been scaled to have a magnitude five times greater relative 
to the scaling used in Figure 6-a (to obtain scaling equivalent to 6-a, reduce 6-b 
by a factor of 5). As the period is increased, we will see a reduction in the fre- 


quency domain's peak energy. Also, as can be seen by comparing Figure 6-b with 
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Figure 5-a, the increased period results in a greater spectral density - there are 


^ more frequency components present. 


Before we increase the period of the pulse train again, look at Figure 6-b and 


note the trend toward grouping of the spectral lines. Also, note the relative ampli- 


. tudes of the grouped spectral lines. Mentally draw a line connecting the peaks of 


the spectral lines and retain this image as we proceed to increase the period of the 
repetitive square pulse. 

^ Let us be really drastic and increase the pulse train period to infinity. When 
we do this we will have a single square pulse as shown in Figure 7-a. The Fourier 
transform magnitude for this single pulse, or if you prefer a pulse train with 
infinite period, is shown in Figure 7-b. How does Figure 7-b compare with the mental 


picture that you constructed from Figure 6-b? 
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a.) A single square pulse. Ь.) Тһе Fourier transform 
magnitude of the square 
pulse in a. 


Figure 7 


What has happened in Figure 7-b is that so many frequency components are 
represented that we no longer see individual spectral lines - the spacing between 
the spectral lines has gone to zero. We see an envelope of frequency content which 
indicates that a single pulse consists of all frequencies in varying amplitude. We 
no longer have the discrete-frequency spectrum associated with periodic waveforms 
but instead have a "continuous spectrum." 

This last concept, the idea of a continuous spectrum, will play an important 
role in the next subject of discussion, the discrete Fourier transform (DFT) and 


its implementation through the fast Fourier transform (FFT). 





The Discrete Fourier Transform 

The Fourier transform that has been discussed up to this point is frequently 
referred to as the continuous Fourier transform (CFT) and is in the form of an 
integral. It is referred to as continuous because the mathematics of the transform 
deal with continuous waveforms. To apply the continuous Fourier transform to a wave- 
form, we must first be able to express that waveform as a continuous function of 
time. For example: a sine wave may be expressed as a continuous function of time 
by f(t) = Asin(2mft). For many of the classic waveforms - sine waves, square 
pulses, etc. - this is a relatively easy matter; however real-life waveforms do not 
generally lend themselves to easy formulation. | 

Fortunately digital techniques, digital computers, and the discrete Fourier 
transform (DFT), which takes the form of a summation, can save us from the mathe- 
matical agony of trying to write and solve formulas for real-life waveforms. The 
discrete Fourier transform is the digital equivalent of the continuous Fourier 
transform and may be used in conjunction with a computer to calculate, from discrete 
samples of the continuous waveform, discrete samples of the continuous spectrum. 

There are two problems related to the implementation of the discrete Fourier 
transform. First, the discrete Fourier transform operates on a sampled version of 
the original waveform and in so doing is subject to errors inherent in the sampling 
technique used. Secondly, and most importantly, computer solution of the discrete 
Fourier transform for any appreciable number of samples is costly in computer time. 
This second item, coupled with the usefulness of Fourier techniques in many appli- 


cations, was the impetus for developing the fast Fourier transform (FFT). 


The Fast Fourier Transform 

The fast Fourier transform is a computer algorithm that allows time-economical 
calculation of the discrete Fourier transform. The time economies associated with 
the FFT arise not out of approximations but rather from a more elegant computational 
approach to the discrete Fourier transform. As compared to straight forward evalu- 
ation of the discrete Fourier transform, the Cooley-Tukey FFT algorithm takes advan- 
tage of a not so obvious approach to calculating the Fourier coefficients with a 
resulting drastic reduction in the number of arithmetic operations required. The 
results obtained by either approach are the same; the FFT simply provides the results 
faster and more economically. 

We should realize by now that the FFT is a fast way to evaluate the discrete 


Fourier transform and that the discrete Fourier transform is a digital version of 
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the continuous Fourier transform. In discussing the continuous Fourier transform, 
several waveform examples were used in exploring the frequency domain. Let us take 
one of these examples now, specifically the continuous sine wave of Figure 2, and 
re-explore the frequency domain through the eyes of the FFT. 

The first thing that we have to understand about the FFT is that it can not 
operate on a continuous waveform. One sine wave example, shown in Figure 8-a, must 
be converted from a continuous function to a digital representation of its former 
self. The Digital Processing Oscilloscope (DPO) is the ideal instrument for doing 
this. Also, the inclusion of the FFT and its inverse in the DPO's software will 


allow us to rely upon this instrument for further discussion of the FFT. 


b.) 
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Figure 8 


In Figure 8-a we have a continuous periodic function of time, a sine wave. We 
are interested in transforming this sine wave into the frequency domain by using the 
FFT. First, we have to obtain and digitize a sample of this waveform. The sample 
should contain enough information, several cycles will be sufficient, to characterize 
the waveform. This requirement is nicely met by the sampling window of Figure 8-b. 

The sampling window of Figure 8-b is the fixed time frame through which we look 
at the continuous wave. When we look through a window at the sine wave, we see only 


the windowed segment and thus have no knowledge of the waveform before or after the 
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window. In effect, we may think of this window as a unity amplitude pulse that is 
multiplied with the sine wave of Figure 8-a. The result of the time domain multi- 
plication, shown in Figure 8-c, is a windowed segment of the continuous sine wave. 
As far as we are concerned, we can say that the sine wave does not exist beyond the ` 
edges of our window. The waveform is no longer periodic. We simply have a small 
segment of a sine wave - a few short wiggles in the infinite continuum of time. 

Now all we have to do is digitize what we see through the window and let a 
computer and the FFT algorithm go to work. When a waveform is stored with the DPO, 
the windowed segment represented by the display is digitized and placed in memory. 
The APD BASIC II software FFT command will cause the minicomputer to transform the 
data to the frequency domain with the results being in complex form, i.e. a real 
part and an imaginary part. To observe the frequency domain magnitude, we can take 
the real and imaginary parts and calculate, with the minicomputer, the square root 
of the sum of the squares. If we have used the "three argument FFT" command, our 
frequency domain magnitude display will be of the positive frequency domain only; 
but remember, the negative frequency domain is the mirror image of the positive 
frequency domain so, for the purpose of display, the positive frequency domain is 
sufficient. 

Let us look at the magnitude of the FFT results for our "windowed" sine wave. 


These results are shown in Figure 9. 





a.) A sine wave in the time domain. Ь.) The FFT magnitude of the sine wave. 


Figure 9 


-12- 


In Figure 9-a we see the sine wave that has been windowed and digitized. When 
we look at its FFT magnitude, we should expect to see a strong indication of its 
frequency. Indeed this is so, as is shown in Figure 9-b. However, there are some 
things about Figure 9-b that may at first be disturbing. For instance, the FFT seems 
to say that there are numerous frequency components crowded about the primary fre- 
quency of the sine wave. At first this may seem to conflict with our intuitive 
expectations, but remember that we windowed the original sine wave which was equiv- 
alent to multiplying by a square pulse. This window or square pulse then might be 
expected to have some influence on the sine wave's frequency domain representation 
and in fact it does. What we are seeing in Figure 9-b is the frequency domain 
representation of the window convolved with the frequency domain representation of 
the sine wave (convolution is a very useful mathematical operation that will be 
covered at a later date). To get an intuitive feel for what this convolution does, 
consider the frequency domain representations of the window and the sine wave before 


windowing. Both of these representations appear in Figure 10. 
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a.) Frequency domain magnitude 
for a rectangular window. 


b.) Frequency domain magnitude 
for a sine wave. 


c.) Result of convolving a with 
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Figure 10 
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further computer calculations. 

The display problem of spectral information being contained within the first 
division is really not a problem at all. Actually, it is a helpful indication! The 
FFT has associated with it the possibility of aliasing errors. Aliasing occurs when 
the sampling rate for digitizing the continuous waveform is less than twice the 
highest frequency present in the waveform. In other words, a waveform is completely 
specified when we sample at a rate sufficient to provide at_ least two samples of its 
highest frequency component. Sampling at a higher rate simply provides extra infor- 
mation which, in the DPO, results in crowding of spectral components to the left side 
of the display. The presence of this crowding, such that the right hand side of the 
display is barren of spectral components, indicates that aliasing or foldover has 
, not occurred. (For more information on aliasing, see APPENDIX A.) 

Probably the major objection to crowding of spectral information to the left 
is that it does not look right. It violates the classical picture of the frequency 
spectrum, and of course, if you are interested in just looking at the spectrum 
magnitude, you would like to see a display that emphasizes the fundamental and its 
immediate harmonics. The solution to the problem is quite simple and serves to 
point out a useful feature of the DPO - variation of sampling rate through the 
horizontal time-per-division control. 

The solution to spectral crowding to the left is to look at more cycles in the 
time domain or if the event of interest is non-repetitive, a pulse, restrict the 
pulse width to the leftmost horizontal division or less. Why does this work? It 
works because the horizontal time-per-div setting effectively controls the DPO's 
sampling rate. Actually the sampling rate itself does not change, but the number of 
samples per cycle for a fixed signal frequency decreases as you increase the hori- 
zontal time-per-div. The DPO takes 512 samples of the displayed time domain wave- 
form. If the display is of one cycle, the waveform will be sampled 512 times per 
cycle; if the horizontal time-per-div control is adjusted to display a greater 
number of cycles of the same waveform, say 25 cycles, then the samples per cycle 
would be 512/25. The result is reduced sampling of higher frequency components, 
remember it takes at least two samples to specify a component, and limiting of the 
FFT display to the lower frequency region. | 

Limiting of the display to the lower frequency region occurs because the Hertz- 
per-div scale factor for the FFT display is derived by dividing 2.56 (in the 4 argu- ` 
ment FFT 5.12 is used) by the time-per-div Scale factor. The display's leftmost 


graticule line is always zero Hertz so, if the scale factor is 2.56/time, you can 
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see that increasing "time" results in changing the frequency domain scale factor to 
exclude the higher frequency range from the display. 

The waveform of Figure ll-a is essentially the same as that in Figure 9-a except 
that in Figure ll-a the horizontal scale factor has been increased by an order of 
magnitude. The result of this increase, and its resulting decrease in sampling, 
becomes obvious when Figure 9-b and Figure ll-b are compared. Note that the fre- 
quency scale factors reflect the order of magnitude difference in the time scale 


factors. 





a.) A sine wave in the time domain. b.) FFT magnitude of the sine wave. 
Figure 11 
The solution is even more dramatically illustrated when the display of Figure 


ll-a is changed to a square wave of the same frequency. Finding the FFT magnitude 


of such à square wave would result in the display of Figure 12. 
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Figure 12 


In Figure 12, if the time domain display had included only one or two Башы of 
the square wave,.we would have seen a resulting crowding to the left. Specifically, 
if the time domain horizontal had been an order of magnitude less, Figure 12 would 
have been scaled to 25.6 KHz instead of 2.56 KHz, and what presently appears in 
Figure 12 would have been displayed entirely within the first horizontal division. 

In spreading the display out, as has been done in Figure 12, the risk of ali- 
asing error increases. Aliasing is not obvious in Figure 12; however the presence 
of a fairly significant frequency component near the right side of the display is 
enough to raise suspicion of more frequency components existing beyond the high fre- 
quency bound of the display. If frequency components do exist beyond the right side 
of the display, they will be "folded" about the right side of the display to appear 
in the display as low frequency aliases. In Figure 12, the right side of the dis- 
play falls midway between the harmonics; therefore if aliasing occurs, the high 
frequency components will be folded back nearly on top of the low frequency compo- 


.nents. For more explicit examples of aliasing, please refer to APPENDIX A. 


Windowing To Improve FFT Results 


Even with the problem of left-side spectral crowding resolved in Figure 12, 
there are some things that might still be bothersome. For example, there appears to 
be, in Figure 12, considerable noise between the harmonic frequencies, not to mention 


the annoyance created by leakage into the side lobes or skirts. 
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The noise problem can be improved upon by signal averaging the acquired signal 
before performing the FFT. The amount of signal averaging required to sufficiently 
reduce the effects of noise depends upon the amount of noise present in the original 
signal. If the original signal is nearly buried in noise, you must signal average 
several thousand times. For signals with moderate noise, as is the case for most 
measurements, 128 averages generally proves sufficient. With low noise, 32 will be 
all that is required. A short APD BASIC program to perform averaging of a waveform 


being stored in DPO location PA could be written as follows: 


100 AVE PA,PB,32 
110 LET PA = PB-MEA(PB) 
120 FFT PA,B,C 

- 130 LET PD - SQR(B*B4C*C) 


Line 100 averages the signal 32 times to remove random noise. Line 110 is a suggested 
step for removing any dc component from the acquired waveform to help expand the dis- 
play vertically. From previous discussion, we are already familiar with the functions 
of lines 120 and 130; they are the FFT and magnitude calculations. The results of 
this program are shown in Figure 13-a and -b. 

. Figure l3-a shows an acquired square wave that has been averaged to remove noise. 
The FFT magnitude of this averaged waveform appears in Figure 13-b. Note that the 
. spectrum display is virtually free of inter-harmonic noise. The problem of noise 
can easily be relieved by signal averaging; however, as can be seen in Figure 13-b, 
the spectrum still includes the skirts, in the display we see primarily an envelope 
formed by the side lobe peaks, associated with the DPO's inherently rectangular 
window. 

The skirts or side lobes associated with the rectangular window arise from the 

sharp discontinuities created by the left and right sides of the window. Note par- 
ticularly the right side of Figure 13-a where the square wave has been abruptly ter- 


minated in mid-cycle by the right side of the windov. 
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а.) Rectangular windowing | b.) Spectrum of the rectangularly 
of a square wave. windowed square wave. 





c.) Triangular windowing d.) Spectrum of the triangularly 
of a square wave. windowed square wave. 


Figure 13. A comparison of rectangular and triangular windowing. 


If the skirts in Figure 13-b arise from the abrupt nature of the rectangular 
window, it is reasonable to assume that a differently shaped window, one with less 
abrupt transitions, would reduce the prominence of the skirts. This assumption is 
born out by Figures 13-c and -d. In Figure 13-c the acquired waveform has been 


multiplied by a triangular window. An APD BASIC program for triangular windowing 
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such as seen in Figure 13-c could be written as follows: 


600 LET EM - 1 

605 FOR J = 0 to 255 

610 LET A(J) = A(J)*J/255 

620 NEXT J 

630 FOR J = 256 TO 511 | 
640 ТЕТ А(Ј) = А(Ј)* (511-7) /255 
650 МЕХТ Ј 

660 ТЕТ ЕМ = 0 


The FFT magnitude of the triangularly windowed waveform, shown in Figure 13-4, 
exhibits skirts that are markedly reduced from those appearing in Figure 13-b. 

A closer comparison of Figure 13-b and -d reveals some differences other than 
reduced skirt prominence. Comparing scale factors and associated spectra reveals a 
reduction of spectrum power for application of the triangular window. This is easily 
understood when one considers that a unity amplitude triangle encloses less area than 
a unity amplitude rectangle of the same horizontal dimensions. In fact, a little 
thought will show that the triangular window contains one-half the area of the rec- 
tangular window, and this difference is seen to be closely reflected in the relative 
amplitudes of the major lobes in Figure 13-b and -d. In Figure 13-d the amplitude | 
is very nearly one-half of that seen in Figure 13-b. 

One other important consideration associated with windowing is illustrated in 
Figure 13-d. The skirts have been reduced by triangular windowing, but what hap- 
pened to the power associated with the skirts? We can not just throw the power 
away, so we must have put it somewhere. A comparative inspection of the major lobe 
widths in Figures 13-b and -d reveals that the spectrum associated with triangular 
windowing contains wider major lobes. Triangular windowing has reduced the power 
in the skirts and redistributed that power to the major lobes with a resulting 
increase in their width. | 

The preceding discussion of windowing should suggest some applications. For 
example, when the situation arises where a low level component resides within and is 
hidden by the skirts of an adjacent high level component, windowing techniques may 
be used to reduce the skirts and reveal the presence of the low level component. 
However, care should be exercised because the possibility exists that the increased 
major lobe width associated with windowing could cause closely adjacent components 


‘to blend together and appear as one. 
Our discussion to this point has dealt only with rectangular and triangular 
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windows. There are many more windowing functions, also called weighting functions 
and convolution kernals, available for preconditioning acquired signals. Table I 
contains some of the more common windowing functions and some of their more impor- 


tant characteristics. 


















Table I 
HIGHEST 
WINDOWING FUNCTION BANDWIDTH SIDE LOBE 
(3 dB) (dB) 
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In Table I the term 8 is used to normalize the theoretical values for general 
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application. Beta, 8, is the reciprocal of the window length in time. For most 
applications concerning the Digital Processing Oscilloscope, beta will be equal to 


the reciprocal of ten times the displayed horizontal scale factor. 


29 


The 3dB BANDWIDTH information in Table I refers to the frequency domain band- 
width of the windowing function's major lobe. The HIGHEST SIDE LOBE (dB) information 
is the amplitude of the highest side lobe in decibels referenced to the peak of the 
major lobe. In looking at these two columns of information, note that, in general, 
as the side lobe level is decreased the major lobe bandwidth increases. In other 
words, the resolution, or ability to distinguish between equal amplitude adjacent 
frequencies, decreases as bandwidth increases. On the other hand, as the side lobe 
energy is decreased the selectivity - ability to resolve adjacent frequencies of 
unequal amplitude - is increased. 

To gain some idea of the relative effects of various windowing functions, you ' 
are referred to APPENDIX B. Here you will find that the same waveform has been 
windowed by each of the functions in Table I and then transformed to the frequency 
domain to illustrate the effects of windowing. Additionally, the APD BASIC programs 
for generating each windowing function are included in APPENDIX B. 

When should windowing be used or when should it not be used? If windowing is 
needed, which windowing function should be used? The answers to these questions are 
summed up by saying, "It depends upon what you are trying to do." If you have adja- 
cent frequency components of nearly equal amplitude, you may not want to window 
because the increasing lobe width may cause them to appear as one component.  How- 
ever if you have a small component near a large one, windowing will reduce the 
skirts of the large one and make the small component visible. The type of window 


used is a compromise between the amount of skirt reduction needed and the allowable 





increase in lobe width. 

The use of windowing and the choice of windowing functions require some prior 
knowledge of the signal to be windowed and some knowledge of what you are looking 
for in the frequency domain. Discussion of windowing to any greater depth would 
exceed our intended scope so we will leave the subject with the following thought. 
Windowing may be thought of as being analogous to the use of light filters in photog- 
raphy. There are some clear cut cases where filters will improve the picture greatly, 
and there are some cases that require experimentation before the fine details are 
arrived at. 

By now you should have a good feel for the Fourier transform. You should be 
aware that it is a means for transforming waveform information from the time domain 
to the frequency domain. You should also be aware that the fast Fourier transform 
operates on a representation of the actual waveform and provides a related frequency 


domain representation; and also that, the act of using a time domain limited or 
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_ windowed representation, along with digitizing, introduces deviation from the ideal. 
In general, these departures create no significant analytic problems, and where 


they do, special techniques usually solve the problem. 
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APPENDIX A 
Aliasing 


Throughout this appendix, a square wave will be used in examples supporting the 
discussion of aliasing. The exclusive use of the square wave is not intended to 
imply restriction of aliasing to this class of waveforms, and it should be realized 
that aliasing is an inherent error associated with sampling in general. 

Before the Fast Fourier Transform can be applied, the analog waveform must be 
digitized, and herein lies the source of aliasing error. To understand what causes 
aliasing, let us begin with the analog waveform shown in Figure A-1 and investigate 


some of its properties and what occurs when they are digitally characterized. 





Figure A-l1. A square wave of peak amplitude A. 


The frequency content of the square wave in Figure A-1 can be examined by 
writing its Fourier series. Any textbook containing mention of Fourier techniques 
invariably uses a square wave as an example for constructing the Fourier series; 
therefore, to expedite the discussion of aliasing, we refer you to any of the many 
textbooks for the details and simply present the results here. The Fourier series 


for the square wave in Figure A-1 is written as follows: 
f(t) » AA/T(sinugt-tl/3sin3ugt-*l/5sin5ugttl/7sin70gt...) Eq-A 


With regard to Eq-A, we should first note that, although only four terms have been 


shown here, the Fourier series is an infinite series. Secondly, we should note that, 
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for the square wave in Figure A-l, the infinite Fourier series consists of sine terms 
only. These sine terms are related such that the square wave is made up of odd har- 
monics with amplitude decreasing as harmonic number increases. (Remember that we 
are discussing a specific example, here, the square wave, and that in general the 
Fourier series is infinite and may contain sine and cosine terms individually or in 
combination.) 

Figure A-2 is a graphical representation of the time domain phase and amplitude 
relationships contained within the truncated series of Eq-A. If we sum the sine 


waves of Figure A-2, the result will appear as shown in Figure A-3. 
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Figure A-2. The first four terms in Eq-A. 
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Figure A-3. The sum of the first four terms in Eq-A. 


'As can be seen by Figure A-3, summing the first four terms of Eq-A has resulted 
in an approximation of the original square wave. The rise time, corners, and flat- 
ness of the reconstituted square wave suffer from the lack of higher frequency compo- 
nents. To improve on this, we can carry the series of Eq-A out to a few more terms 
and add these to those already incorporated in Figure A-3. Doing so would result in 
a waveform with faster rise time and reduced ripple - a closer approximation to the 
original square wave of Figure A-l. If we continue to do this, add more terms, we 
will finally reach a point, far short of having an infinite frequency term, that 
produces a very close approximation to the ideal square wave. Such an approximation 
is generally referred to as a real-life square wave, and we can say that the remaining 
terms necessary to the ideal square wave have been attenuated to a level too low for 
consideration or have gone to zero. 

It can be seen from the preceding discussion that a real-life square wave can 
essentially be described by a Fourier series wherein the higher frequency terms have 
been attenuated to a level of zero. Another way of looking at this is to consider 


the hypothetical situation where we generate an ideal square wave. As soon as we 
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try to transmit this square wave to an amplifier or some other instrumentation, 

we find that the square wave suffers some degree of degradation. This will occur 
because even the best transmission mediums and instrumentation have bandwidth limi- 
tations. Frequency components that exceed the bandwidth of the system will be 
significantly attenuated. In short, life is a low-pass filter. 

Realizing now that a nonsinusoidal waveform is in essence the sum of a finite 
number of sinusoidal terms, we can get on with the task of digitally characterizing 
the waveform. Digitizing an analog waveform requires that we sample the waveform's 
amplitude at enough time intervals to define the waveform. The number of times that 
any waveform is sampled in a fixed time period is referred to as the sampling rate, 
and the Nyquist criterion for completely defining a waveform requires that the 
sampling rate be twice the highest frequency present in the waveform. For example, 
if the highest frequency term in a waveform is known to be the thirtieth harmonic, 
a sampling rate of twice the frequency of the thirtieth harmonic will be sufficient 
to completely define the waveform. If the sampling rate is reduced to a lower fre- 
quency, aliasing or foldover will occur. 

As a matter of definition, the frequency equal to one half of the sampling 
frequency is called the Nyquist or folding frequency. If the sampled waveform con- 
tains frequency components above the folding frequency, those high frequency compo- 
nents will impersonate, or assume as an alias, a frequency below the folding fre- 


quency. To visualize the mechanics of this impersonation, refer to Figure A-4. 


Figure A-4. Insufficient sampling of a high frequency component will result in a 
low frequency alias. 
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To set the stage in Figure A-4, let us assume that the ten cycles of sine wave 
shown with the solid line represent a high frequency component, say 100 KHz, of a 
waveform that is being digitized. Let us further assume that the heavy dots repre- 
sent digital amplitude samples relative to the analog-to-digital reference, shown as 
a solid horizontal line. Note that, for the example in Figure A-4, there are twelve 
samples for ten cycles of sine wave. If the sine wave has a frequency of 100 KHz and 
is sampled 1.2 times per cycle, it follows that the sampling frequency is 120 KHz. 
We know from past discussion that this sampling rate is too low for a 100 KHz compo- 
nent, the aliasing or foldover frequency is 120 KHz/2 or 60 KHz, and as a result we 
should see aliasing. 

The low frequency that is impersonated by the insufficient sampling of the high 
frequency component is shown in Figure A-4 as the dotted sine wave. Note that the 
impersonation passes through each of the sample points. A little further inspection 
shows that, for the assumed sampling frequency, the dotted sine wave has a frequency 
of 20 KHz. 

If we re-examine the frequencies discussed relative to the example in Figure 
A-4, we can begin to see a pattern in aliasing. We assumed a high frequency compo- 
nent of 100 KHz. From the arrangement of digital samples on this high frequency 


component, we were able to determine that the sampling frequency was 120 KHz with a 


corresponding folding frequency of 60 KHz. For this case, the folding frequency is 
40 KHz below the high frequency component. We also noted in Figure A-4 that the 


sampling points also describe points on a low frequency sine wave. This low fre- 
quency sine wave, the alias, was determined to have a frequency of 20 KHz - it occurs 
at a frequency 40 KHz below the folding frequency. Do you see the pattern? 

When a sampled waveform has frequency components above the Nyquist or folding 
frequency, those components will be folded about the Nyquist frequency into the fre- 
quency domain below the Nyquist frequency. To aid in visualizing this folding action, 


refer to Figure A-5. 
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Figure A-5 


In Figure A-5, the frequency domain magnitude in solid lines represents those 
portions below the folding frequency. The dotted portion represents the spectrum 
above the folding frequency. Foldover appropriately derives its name from the 
folding of the high frequency spectrum into the low frequency spectrum. When the 
spectrum above the folding point contains significant components, those components 
will appear as aliases below the folding point. 

Aliasing can be demonstrated with the Digital Processing Oscilloscope by gener- 
ating a near ideal square wave and then obtaining its FFT magnitude. Square waves 


can quickly be generated by typing in the following: 
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LET A = .6283 

INT A,PA 

LET PB = SIN(N*A) 

LET C = SGN(B):LET C(O) = 1 
LET PC = C 


where N is set to the number of cycles desired. Once the basic square wave is 
generated, the vertical and horizontal scale factors and units may be set as desired. 
The FFT magnitude of 24 displayed square-wave cycles, with cosine windowing to reduce 
skirts, is shown in Figure A-6. 
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Figure A-6. The circled components are aliases. 





In Figure A-6, the folding frequency is the right side of the display and the 
aliases are circled for easy identification. We can in general determine the pres- 
ence of aliasing and identify the aliased frequencies by a few rules of thumb. 

We know from previous discussion that a square wave is made up of harmonics 
that are integer multiples of the fundamental and decrease in magnitude as frequency 
increases. By observing the first three uncircled somponeute in Figure A-6, we can 
determine the correct spacing for the harmonics. With this frequency increment or 


spacing in mind, we can predict each successive harmonic's location as we travel to 


\ 
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the right in Figure A-6. When we reach the last uncircled component to the right, 
the predicted location of the next component will fall beyond the edge of the dis- 
play. We know from the general trend of magnitudes that this component will have 
a decreased but still significant magnitude. We also know from its predicted loca- 
tion that the component exists beyond the right side of the display by "x" distance; 
hence it should be folded about the right side of the display to appear "x" distance 
to the left. This component does appear at the predicted location and is shown as 
the circled component closest to the right side of the display. If we continue to 
the left from this first alias with the same spacing predictions, we will find the 
locations of succeeding aliases; and if after we reach the left side of the display 
and we predict further significant components, we will see these components folded 
about the left side of the display. In Figure A-6, the components enclosed by small 
circles were folded about the left side of the display. They are aliases whose real 
frequencies are greater than twice the folding frequency. 

In summary, we can state some basic guidelines for identifying aliasing in the 


Digital Processing Oscilloscope. 


1) If the FFT magnitude display has significant components near 
the right side of the display or does not go to zero and remain 
at zero thereafter, aliasing has probably occurred. 


2) If there are equally spaced components present in the FFT 
magnitude that have decreasing magnitude from right to left 
(the large circles in Figure A-6), then these components are 
probably aliases. 


3) If aliases can not be identified and yet there are significant 
components occurring near the right side of the display, check 
the spacing between the lower frequency components and compare 
this to the distance from the right edge of the display to the 
rightmost displayed component. If the ratio of the two distances 
is two to one, it is likely that the high frequency components 
are folding back on top of the low frequency components. 

(Figure 12 in the text is an example of where high frequency 
components have folded back almost on top of the low frequency 
components.) 


4) Low level spectral components residing between higher level 
spectral components are not always aliases. If these components 
are located exactly half way between the adjacent higher level 
components, they may be valid harmonics. In the case of a non- 
symetric square wave, the degree of nonsymetry is indicated by 
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the magnitude of the even harmonics. In a nearly symetric 
Square wave the magnitude of the even harmonics will be 
very small relative to the magnitudes of the adjacent odd 
harmonics. 


As was stated at the beginning, these are just guide lines and should not be 
considered all encompassing or binding. In fact, these guide lines are primarily 
restricted to periodic waveforms. Aliasing in the case of nonperiodic waveforms is 
more difficult to detect; however aliasing may be assumed to exist if the frequency 
domain magnitude of the nonperiodic waveform has not gone to zero before the right 
side of the display. 

Since aliasing arises from insufficient sampling of the original waveform, it 
, is obvious that the cure for aliasing problems is to assure sufficient sampling of 
the original waveform. In many cases the cure is to simply digitize fewer cycles 
of the time domain waveform. Unfortunately, there are a few cases where the orig- 
inal waveform is not high frequency limited (i.e. very fast rise times, responses 
from high-pass filters, etc.). In these cases, aliasing can be prevented by fil- 
tering before digitizing. Filters used for this are referred to as antialiasing 
filters and are designed to limit the high frequency content of the filtered wave- 
form to a known and acceptable cutoff frequency. 

For most Digital Processing Oscilloscope applications, aliasing is not a 
problem. For Hose Gases where aliasing does occur, appropriate measurement 


technique can usually eliminate or minimize the problen. 
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APPENDIX B 
Windowing 


Contained within this appendix are several examples of windowing. These examples 
are arrived at by always windowing the same waveform so that the relative effects of 
each window may be seen in the FFT magnitude display. Each of the eight windowing 
examples consist of a display of the windowed waveform, a display of the FFT magni- 
tude, and the APD BASIC II, V01-03, program for windowing a waveform stored in the 
TEKTRONIXe P7001 Processor memory location A. This appendix also contains a pro- 
gram for demonstrating windowing and a brief discussion of the program call button 
assignments. 

The waveform used was derived by feeding the DPO's 4V CALIBRATOR output into 
the low-pass filter of Figure B-1. The output of the low-pass filter was signal 
averaged and then a sine wave was added to the waveform. This sine wave was of fre- 
quency slightly greater than that of the square wave and had an amplitude equal to 
- 0.09 of the square wave's amplitude. The low-pass filter was used to prevent ali- 
asing, the signal averaging was done to remove noise, and the sine wave was added 
to aid in comparing the effects of various windows. 

The most significant items to note in the FFT magnitudes shown are: the pres- 
-ence of a low-level component that is closely adjacent to the fundamental and the 
effects of various windows on this component, the relative levels of each FFT magni- 
tude for various windows, the relative effect of various windows on the skirts or 


side lobes, and the relative major lobe widths for each window. 
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Windowing Demonstration Program 


AVERAGE FA. PE. 122 

(ЕТ PA=E-MEACE > 

(ЕТ РЕ-РН 

STOP 

LET FA=FPB: HOLD. A 
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Program Call Button Assignments 


Program 
Call Button Function 
1 Averages PA into PB 128 times, removes the 
mean, and stores the waveform in both PA 
and PB. 
2 Triangular Places the unwindowed 
3 Cosine or original waveform 
in PA, generates and 
À Gogine squared places the selected 
5 Cosine cubed windowing function in 
6 Cosine to the fourth Poa and wiadows the 
original waveform by 
7 Hamming the selected windowing | 
8 Triplet function. 
| 
9 Places the contents of PB, the original or | 
unwindowed waveform, in location PA. n 
10 | Calculates the FFT magnitude of the wave- | 
form in PA and places the result in PB. | 
11 Generates a sine wave and adds it to the | 


waveform in location PA. The result is 
stored in both PA and PB. 
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Initiating And Running The Program 
Set the time base to 1 mS/div and vertical amplifier to .1 V/div. 


Feed the 4 V CALIBRATOR signal through the low pass filter shown in Figure В-1. 
Store the filtered signal in location A. You should see ten cycles of the 
filtered CALIBRATOR signal. The filtering is necessary to prevent aliasing.* 


Press program call button number 1 to signal average the acquired waveforn. 


Press program call button number 11 to add a sinusoidal component that will be 
close to the fundamental of the filtered CALIBRATOR signal. If you wish to vary 
the added component's proximity to the fundamental, change the value of N in 
line 1105 of the program. To change the component's amplitude, change the value 


of the constant in line 1115 of the program. 


Press program call button number 10 to obtain the FFT magnitude. Do you see the 


added component near the fundamental? 


Press program call button number 2 for triangular windowing and button 10 for 


the FFT magnitude. Can you see the added component? 


For other windows, press program call buttons 3 through 8. After selecting a 


window, press button 10 to obtain the FFT magnitude. 


Program call button 9 is used to place the unwindowed waveform back into PA. 


*For the purpose of display, ten cycles of the CALIBRATOR waveform are acquired. 

This provides an FFT magnitude display that is spread out for visual observation; 
however the acquisition of ten cycles of CALIBRATOR waveform results in undersampling 
and thus aliasing. Though not apparent in the FFT of the rectangularly windowed 
waveform (rectangular windowing is inherent in the DPO), other windows tend to accen- 
tuate the high frequency aliases through side lobe reduction. To maintain an 
uncrowded display without the clutter of high frequency aliases, low-pass filtering 
is necessary. 
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